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Abstract 

A discrete iV— point Runge-Kutta version H^ N \X) of one of the simplest 
non-Hermitian square-well Hamiltonians with real spectrum is studied. A 
complete set of its possible hermit izat ions (i.e., of the eligible metrics Q^ N \X) 
defining its non-equivalent physical Hilbert spaces of states) is constructed, 
in closed form, for any coupling A G (—1, 1) and any matrix dimension N. 



1 Introduction 



1.1 Bound states in Runge-Kutta approximation 



The concept of the solvability of a dynamical model in physics is rather vague. 
Its definition is usually adapted to the range of expected applications. By 
some authors even the single-particle motion along a finite one-dimensional 
interval would be called solvable only if the underlying ordinary differential 
Schrodinger equation for bound states 



dx 2 



ip{x) + V(x) ip(x) = E ip{x) 



ip(±L) = 



(1) 



proved reducible to the Gauss' or confluent hypergeometric equation. Re- 
markably enough, even this extremely narrow specification of solvability finds 
very plausible physical reasons in the related shape invariance of potentials 
V(x) and/or in their close relationship to supersymmetry pQ. 

In what follows we shall rather accept different terminology in which one 
treats virtually any sufficiently smooth and real potential V(x) in Eq. (Tj[|) 
as exactly solvable, i.e., solvable, in principle, with arbitrary precision, in 
a purely numerical setting at least. Typically, the Runge-Kutta (RK, [2]) 
discrete approximation of the above equation, viz., 



i>{x k -x) - 2 ip(x k ) + ip(x kA 



h 2 

0,±1, 



+ V(x k ) ip{x k ) = E ip{x k ) 



±K 



±L. 
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may be used to reduce the original bound-state problem to a routine computer- 
assisted diagonalization of the finite-dimensional RK matrix Hamiltonian 
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\ 



2 + h 2 V{x_ 1 ) -1 

-1 2 + h 2 V(x ) -1 

-1 2 + h 2 V(x 1 ) 



\ 



(3) 



/ 



The choice of the lattice distance h 1 is only dictated by the required 
precision of reproduction of the original energies and/or wave functions. 



1.2 VT— symmetric Runge-Kutta models 

The standard formalism of quantum theory admits Schrodinger Eqs. (CQ) 
and/or (T5]) which generate the real bound-state spectra from certain complex 
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potentials V(x). Conventionally, these potentials are called VT — symmetric 
(cf. reviews [3], H] or Appendix A for more details). For a sensible extension 
of the concept of solvability to this new context it is most important that the 
VT— symmetric models are characterized by the Hamiltonian-dependence of 
their physical Hilbert spaces 7i^ s \ Although the superscript ^ stands here 
for "standard" , the inner product is defined in a nonstandard manner in this 
space, by a formula containing an ad hoc metric operator O = Q(H) ^ I (cf. 
Eq. ( 1271) in Appendix A below). 

In spite of the existence of several powerful techniques of reconstruction 
of O for differential Eq. (JT|) [3], [5] one can rarely find a closed-form result (for 
illustration check a few samples in Refs. [Hl[7]). Many successful constructions 
rely upon various assumptions requiring, e.g., the existence of a charge of the 
system [HIE] or of some other additional and/or complementary observable(s) 
O [TU]. Moreover, people are mostly able to obtain 6 just in an approximate 
form, say, of perturbation series [TT] . 

Some of these difficulties have been addressed in our papers [12]. We 
restricted our attention to the discrete Schrodinger Eq. ([2]) considered at 
a pre-determined precision, i.e., at a fixed spacing constant h > 0. This 
enabled us to broaden the scope of the theory and to consider certain nonlocal 
generalizations of interaction terms. More explicitly, we complemented the 
diagonal elements V(xk) in the RK matrix Hamiltonian ([3]) by a set of real 
chain-coupling constants in a way which still left the resulting asymmetric 
real matrix H tridiagonal, 



/ 



H 



\ 



2 + h 2 V(x^] 
-1 + u 



-1 - u 
2 + h 2 V(x ) 

-1+Ux 



-1 — U\ 
h 2 V{x X ) 



\ 



(4) 



For several special cases of this family of Hamiltonians we then constructed 
particular diagonal-matrix metrics O ^ I in closed form. 

This partial success of the project encouraged us to re-open the question 
of a complete solvability of a VT — symmetric model in our subsequent paper 
|13j . In the cryptohermitian scenario characterized by the nontriviality of 
the metrics 6 ^ I we found a complete set of metrics for a special case 
of Eq. (J4]) with V(xk) = (i.e., without any complex local force) and with 
U\ = u 2 = ■ ■ ■ = and u_i = m„ 2 = . . . = 0, i.e., with the single off-diagonal 
real coupling constant g = Uq ^ representing a nonlocal, manifestly non- 
Hermitian potential. 

In spite of the feasibility of such a construction we still felt disappointed 
not only by the necessity of the really lengthy calculations but also by the 
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comparatively complicated structure of matrix elements of metrics G. Al- 
though these elements were expressible in terms of closed-form polynomials 
in g, the degree of these polynomials grew quickly with the cut-off dimen- 
sion N of the g— dependent Hilbert space 7i^ s \ In this sense the results 
of Ref. [13] proved discouraging, especially in the light of a really extreme 
simplicity of the underlying single-center interaction. 

In our present paper we intend to report a return to optimism. Firstly, in 
a preparatory Sec. [2] we shall introduce a double-center model and show that 
its more complicated dynamics does not worsen the feasibility of calculations 
nor a guarantee of the reality of the energy spectra. Next, we shall formulate 
our project of construction of all the eligible metrics for this model in Sec. [31 
In contrast to the similar results of Ref. [13J we shall be able to show here 
that the present version of the VT— symmetric discrete square-well model is 
much more friendly since the related menu of metrics B exhibits a paradoxical 
decrease of complexity during the increase of dimension N. In Sec. H] (dealing 
with exceptional, one- or two-diagonal metrics) and Sec. [5] (describing all the 
remaining metrics) this fact will enable us to find and prove results valid at 
all N. A compact rigorous proof will be delivered confirming the validity of 
our closed and elementary explicit formulae for metrics 0^(A) at any given 
dimension N and coupling A G (—1,1). 

2 Discrete non- Hermit ian square wells 

The key motivation of the present paper resulted from our study of Refs. [H] 
and p] where the differential version ([TJ) of a square-well Schrodinger equa- 
tion has been analyzed with a VT — symmetric point interaction localized in 
the origin (i.e., at x — 0) and at the two distant points (x = ±L), respec- 
tively. The former, simpler, single-center arrangement found its discrete RK 
analogue in the models of Ref. [T3] . Our present paper will offer the discrete 
RK complement to the latter, more ambitious study [7|. 

2.1 Hamiltonians 

We shall analyze the two-center boundary-interaction N— dimensional mod- 
els H( N \\) which form the family 
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(etc) and which have the following tridiagonal N by N matrix form in general, 
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We intend to pay attention to these models in bound-state regime. Hence, 
we have to clarify, first of all, the structure of the domain of couplings A for 
which the energy spectrum remains real. 

2.2 The reality of the spectra of energies 




Figure 1: Real spectrum and its degeneracy and complexification at N = 4. 

In the first step, at N = 3 we obtain the three easily evaluated eigenvalues 
E = 2, E ±1 = 2 + (2-2A 2 ) 1/2 

5 



which are all real and non-degenerate inside the interval of A G (—1, 1). Next, 
at iV = 4 we arrive at the four eigenvalues 

E ±l/2 = 3/2 ± 1/2 (5-4 A 2 ) 1 / 2 , E ±3/2 = 5/2 ± 1/2 (5-4 A 2 ) 1 / 2 

which stay real within a larger interval of A G (— v5/2, y/E/2). This may 
seem to indicate that the domain of the admissible values of A may vary with 
the dimension. Fortunately, it is not so. Figured] clarifies the apparent puzzle 
by showing that the condition of non-degeneracy gets violated precisely at 
the boundary of the same, dimension-independent interval A G (—1, 1). 
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Figure 2: The A— dependence of energies at N = 5. 

Proceeding to N = 5 we reveal that the five eigenvalues 

E = 2, E ±1 = 2 ± (1 - A 2 ) 1/2 , E ±2 = 2 ± (3 - A 2 ) 1/2 

are all real and non-degenerate inside the same interval of A G ( — 1,1). Our 
Figure [2] displays the coupling-dependence of these energies in full detail. 

In the next case, at N = 6, the closed formulae for the energies become 
clumsy. Still, Figure [3] demonstrates clearly that the six eigenvalues behave 
in expected manned at all A G (—1, 1). At iV = 7 (cf. Figure H]) the formulae 
become, paradoxically, simpler, defining all the spectrum by the compact 
equations 



E = 2, E ±1 = 2±l/2y / 8-2A 2 -2v / 8 + A 4 , E± 2 = 2 ± V2 - A 2 , 

E± 3 = 2 ± 1/2 ^8 - 2 A 2 + 2 V8 + A 4 . 
Next, two separate and rather complicated equations of fourth order de- 
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Figure 3: Real parts of energies at N = 6. 




Figure 4: The A— dependence of energies at iV = 7. 



termine the spectrum at even N = 8 (cf. Figure [5]). In contrast, our last 
illustration at iV = 9 (cf. Figure [6]) yields closed formulae again, 



E = 2, E±i = 2±l/2y6-2A 2 -2v / A 4 -2A 2 + 5, 



E± 2 = 2 ± 1/2 \l 10 - 2 A 2 - 2 VA 4 + 2 A 2 + 5 
E± 3 = 2±l/2y6-2A 2 + 2 V\ 4 - 2 A 2 + 5 , 



£ ±4 = 2 ± 1/2 y 10 - 2 A 2 + 2 VA 4 + 2 A 2 + 5 . 

This indicates that the secular polynomials are simpler at odd dimensions. In 
all the Figures [1]- [6] the interval of the allowed couplings A G (—1, 1) does not 
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Figure 5: The A— dependence of energies at N = 8. 

vary with the growth of the dimension N. An independent confirmation of 
this feature of our model will follow, later, from the existence and invertibility 
conditions imposed upon the metric at any N. 
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Figure 6: The A— dependence of energies at N = 9. 

Our pictures illustrate that the extreme values of A = ±1 correspond 
to the triple confluence of the energies E±i with Eq at odd iV and, in our 
notation, to the incidental degeneracy of E + i/ 2 with E_ 3 / 2 at even N. In 
the language of Refs. [T5] one can conclude that at any dimension N > 4 
our model possesses strictly four fragile energies which complexify during 
the transition from the Hermitian limit A = to the asymptotic, strongly 
non-Hermitian regime under very large |A| 3> 1. 
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3 Hermitization 



The study of non-Hermitian Hamiltonians H with real spectra is easier when- 
ever there exists an invertible map Q of H upon an isospectral operator (say, 
upon 1) = QHQ^ 1 ) which is Hermitian (it is, sometimes, called Dyson's 
map). The latter operator f) (acting in some abstract physical Hilbert space 
H.( p \ cf. [UE3IIIIIS]) ma y be assumed prohibitively complicated (otherwise, 
there would exist no good reason for studying H). The space may be 
assumed endowed with the usual Dirac's trivial metric 0^ = I. By defini- 
tion, the latter space must be unitarily equivalent to its equally acceptable 
physical alternative Ti.^ [T7] . 

In applications, the role of the Dyson map Q degenerates to the intro- 
duction of the elementary product = called metric in TC^ S '. For this 
reason we shall only be interested here in the constructions of the matrices 
6 = Q(H). 



3.1 The matrices of metrics 

In the theoretical setting outlined in Appendix A below, the key questions 
to be answered in connection with the analysis of models sampled by Eq. (jSJ) 
result from the ambiguity of the assignment of the metric to a given Hamil- 
tonian H [51 [TU]. For the purposes of clarification of the roots and forms 
of this ambiguity the Runge-Kutta discretized Hamiltonians prove particu- 
larly suitable. Indeed, in the related finite-dimensional (i.e., N— dimensional) 
standard physical Hilbert spaces Ti,^ the admissible metrics = Q(H) do 
only form a finite, strictly N— parametric family. Thus, we may decompose 

N 

ew = X>*f° (6) 

k=l 

and require that the individual Hermitian components of our metric are 
some extremely simple pseudometric (i.e., not necessarily positive definite) 
matrices with, presumably, sparse-matrix structure. 

From the purely algebraic point of view the correspondence between H 
and 0(if) is exclusively specified by Eq. ( |25|) of Appendix A which can be 
rewritten in the explicit linear algebraic form 



N 

E 

fc=i 



[(H^). k e kn -Q jk H kn \ =0, j,n = 1,2,..., N. (7) 



At the first sight, the direct use of such a systems of iV 2 equations for the 
determination of the matrix elements of = Q(H) looks discouragingly dif- 
ficult. Fortunately, not all of these equations are linearly independent. The 
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number of unknowns is also lowered by the necessary Hermiticity of accept- 
able matrices = Q\ Still, for any given Hamiltonian H, an encouragement 
and insight into the generic structure of the solutions = Q(H) may only 
be acquired step by step, by the patient solution of Eq. (j7|) starting from the 
smallest dimensions N. 

Our first nontrivial real and asymmetric, i.e., VT— symmetric and non- 
Hermitian Hamiltonian H^ 3 \\) has already been studied in the different 
context (viz., in connection with the cryptounitary description of scattering, 
cf. section 2.1 of Ref. [IE])- In the present, bound-state version of this 
model we may expect that any eligible H^ 3 \ A) — dependent metric G^ 3 ^(A) 
compatible with Eq. ([7j) will have the real and symmetric six-parametric 
matrix form 

a b c 

e= b f g 

c g m 

The values of its six free real parameters a — m are only restricted by the 
sequence of nine linear relations (jTJ) and by the positivity requirement > 0. 
It is easy to verify that the six nontrivial items of Eq. ([7j) number 2, 3, 4, 6, 7 
and 8 read, respectively, — / + /A + a + aA + c + cA = 0, — g + g\ + b — b\ = 0, 
— a— aA— c— c\+f— f\ = 0, —c—c\—m—m\+f—f\ = 0, —b+b\+g—g\ = 
and — / + f\ + c + c\ + m + m\ = 0. At A G (—1, 1), both the items 
number 3 and 7 give b = g while the comparison of items number 2 and 8 
gives m = a. The remaining four nontrivial equations all coincide with the 
constraint /(l — A) = (a + c)(l + A). We may summarize that the complete 
solution of Eq. has the three-parametric form 



e 



(3) 

(o,6,c) 



(A) 



a b c 

h (a+c)(l+A) , 

u 1-A u 

c b a 



(8) 



At c = we get the tridiagonal matrix 



a 
b 




b 

1-A u 

b a 



which becomes strictly diagonal at b = where one could also set, without 



any loss of generality, a 



a 



(1-A)/(1 + A) yielding 



6 



(3) 

a,0,0) 



(A) 



a 







1 








a 



a 



1 + A 



(9) 



At any nonvanishing A G (—1,1) the latter, diagonal metric is safely positive 
and manifestly non-Dirac, 7^ /. At the same time, the positivity of its non- 
diagonal predecessors must be guaranteed by an explicit (and not too easy) 
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specification of the admissible domain of parameters a, b and c. Without 
this guarantee we may only speak about potentially non-invertible or indef- 
inite pseudometrics replacing, accordingly, also the symbol Oj^ c \(A) for 

the metric, say, by a less specific symbol Q^\ c )(A) whenever appropriate 
(cf., e.g., Appendix B). 



3.2 The set of simplified pseudometrics 

We shall see below that the feasibility of transition to higher dimensions 
N > 3 will be rendered possible by the Runge-Kutta algebraization of the 
Hamiltonian as well as by our localization of interaction far from the origin. 
We shall reveal that these features of our model open the path toward an 
enhancement of efficiency of the construction of via expansion ([6]) where the 
individual pseudometric components V may be sought in certain simplified 
sparse-matrix forms. Naturally, Eq. (j7j) may be required to be satisfied also 
by every individual component matrix V = , 



N 



V H 



1,2, 



N. 



(10) 



The explicit solution of the latter, simplified system of equations will be 
further facilitated by the tridiagonal matrix structure of our Hamiltonians. 
This can be illustrated at N = 4 for which the real, VT— symmetric and 
non-Hermitian Hamiltonian H^>{X) admits the following real and symmetric 
ansatz for the metric 



e (4) (A) = e 



(4) 

(a,b,c,d) 
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b 
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f 
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n 
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J 



Out of the sixteen linear relations in ((Tj) only the ones with numbers 1, 6, 
11 and 16 are trivially satisfied. Further, items 3 and 13 give h = c, items 
7 and 10 yield m — f while the comparison of 2 with 15 gives j = a and 
the comparison of 8 with 9 yields n = b. The rest of the set of constraints 
degenerates to the doublet of requirements / = 6+ ^_ 1 ^ A - > and g 



b+d(l+X) 
l-A 



so that the exhaustive solution of Eq. (j7]) has the four-parametric form as it 
should, 
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(4) 



(A) 
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b 
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d - 
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c+a(l+ 
1-A 


A) 


b+d(l+X) 
1-A 


C 


c 


b+d(l+ 
1-A 


A) 


c+a(l+A) 
1-A 


b 


d 


C 




b 


a 



(12) 
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At d = this matrix becomes pentadiagonal while the additional constraint 
c = makes it tridiagonal. Finally, at b = we get the diagonal metric 



Q (4) 

w O,o,o,o) 



(A) 



a 














1 














1 














a: 



A 



a: 



1 + A 



(13) 



Up to two elements the latter metric coincides with the most common Dirac's 
Q(Dirac) = L one can expect (and we shall verify below) that the existence 
of the similar diagonal metric is a generic feature of our model at all the 
dimensions N = 3, 4, 



4 Exceptional metrics 

The well known ambiguity of the correspondence between Hamiltonian H 
and metric 0(H) [10] seems to be particularly important when one tries to 
improve the precision of RK approximation and to guess the form of metric 
q( n + 1 ) from the knowledge of Q( N \ For such an extrapolation the complete 
solution of Eq. (J7j) at the smallest dimensions A^ < 4 does not suffice. More- 
over, just a marginal help may be extracted from the related mathematical 
literature (dating back to the early sixties [16]). More insight has only been 
obtained during the first use of the metrics 6 ^ I in nuclear physics [9j [10] 
and during the development of VT— symmetric quantum mechanics. In the 
latter setting an additional constraint has been accepted and the metric 
has been assumed factorized into the product of parity (i.e., of one of the 
pseudometrics V) with the so called quasiparity [19] or charge [8]. 

Even under the latter class of additional assumptions the metric 0(H) 
may remain non-unique [5j M, HH EU] . This means that the usual choice of the 
Hamiltonian H should be accompanied by the additional phenomenological 
or pragmatic considerations, i.e., by the physics-dictated or comfort-dictated 
specification of some optimal or exceptional hermitizing metrics 0(H). 

4.1 The strictly diagonal metrics 

In section [3] three important generic features of our present square-well model 
have been revealed. Firstly, the eligible matrices were found to possess 
more symmetries than expected. Secondly, the structure of these matrices 
appeared to simplify when one keeps just single first-line matrix element 
different from zero. Thirdly, the simplest form of matrix elements seems to 
be achieved when the first-line matrix elements a,b, . . . are properly rescaled 
in a coupling-dependent way. 
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In practice, the most prominent role will always be played by the matrices 
9 ^ / which are as close to diagonal ones as possible having a (= O n ) ^ 
while b (= G12) = c (= 613) = . . . = 0. In this sense our first present 
interesting result is the following one. 

Theorem 4.1. For any matrix Hamiltonian (EP with coupling A G (—1,1) 
and dimension N — 3, 4, . . . t/iere always exists a diagonal, positive definite 
metric matrix which differs from the Dirac's Q = I just at the elements Qkk 
with k = 1 and k = N. Its explicit form is 



U (a,0 v ..,0) 



(A) 



a 









1 
















1 










a 



a 



a{\) 



1 - A 
1 + A 



(14) 



Proof. Equation ( jl~4|) coincides with Eq. ([9]) at iV = 3 and with Eq. ( {TBI at 
N = 4. At any integer N > 4 the insertion of ansatz (fT4|) converts Eq. (J7J) 
into identity. □ 

Corollary 4.2. T/ie spectrum of Hamiltonian (TJP coupling A G (—1, 1) 
zs real at any dimension N = 3, 4, . . .. 

Proof. The existence of the metric guarantees that in the corresponding 
finite-dimensional Hilbert space the Hamiltonian matrix is Hermitian. □ 

Just an inessential modification of the above construction leads also to 
the following interesting observation. 

Proposition 4.3. For any matrix Hamiltonian (TJJ) with coupling A G (—1, 1) 
and dimension N = 3,4,... there always exists an antidiagonal, indefinite 
pseudometric matrix of the form 



Q 



(JV) 

(0,0,...,0,a) 



(A) 








a 







1 






1 







a 







a: 



a(A) 



1 - A 
1 + A 



(15) 



Proof. At iV = 3, equation ( |T5l) degenerates to formula (JSj) at a = b = and 
c = a. At N = 4, equation ( fT5l) degenerates to formula ( |T2l) ata = 6 = c = 
and <i = a. At any higher integer N > 4 the insertion of ansatz (|T5|) converts 
Eq. ([7]) into identity. At the same time one can easily verify that matrix (fT5|) 
is not positive definite at any N >2. □ 
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Remark 4.4. It is worth noticing that before one performs the continuous 
h — > limit of the present Runge-Kutta picture, certain metric operators 
obtained in Ref. [5J at A = may be interpreted as positive definite super- 
positions of our two one-diagonal matrices flHJ) and (IT5|) . 



4.2 The strictly bidiagonal pseudometrics 

In our present concrete realization of expansion (jHJ) the general matrix of 
metric 0(1 \(A) will always be written as a linear superposition of the 



'(a,b,-) 

following iV— plet of specific pseudometric matrices 



2(2,.,o)(^) > 



?f°(A)~Q 



(JV) 

(0,6,0,. ..,0) 



(A) 



Their mutual independence is trivially guaranteed by the presence of the 
mere single nonvanishing matrix element in their respective first lines. 

On the level of our above-mentioned N < 4 experience and constructions 
we just managed to guess and prove the generic form of the diagonal and 
antidiagonal pseudometrics Pi (A) [cf. Eq. (1141) ] and 7-V(A) [cf. Eq. (TI5]) ]. 
Further insight must be built using the next, N = 5 ansatz for the real 
and symmetric metric 0®. Using symbolic manipulations we managed to 
reveal that without any loss of generality this ansatz may be written in the 
simplified, symmetrized nine-parametric form 



e 



(5) 

(a,b,c,d,e) 



abode 
b f g h d 
c g m g c 
d h g f b 
e d c b a 



(16) 



Out of the twenty five independent relations (j2j) we managed to eliminate nine 
trivial ones yielding the four quadruplets of equations which enabled us to 
eliminate m = f+h—c(l + X) (from items 8, 12, 14 and 18), g = (b+d) /(l— A) 
(from items 3, 11, 15 and 23), h = (c + e(l + A))/(l - A) (from items 4, 10, 
16 and 22) and / = (c+ a(l + A))/(l - A) from items 2, 6, 20 and 24. Thus, 
we arrived at the complete five-parametric solution @[a\ crfe )(A) of Eq. ([7]), 



d 



b 

c+q(l+A) 

1-A 

b+d 

1-A 
c+e(l+A) 

1-A 

d 



c+q(l+A) 
1-A 



C 
b+d 
1-A 



+ ^^-c(l + A) f±f 



b+d 
1-A 

C 



(a,b,c,d,e) ' 
d 

c+e(l+A) 
1-A 
b+d 
1-A 
c+q(l+A) 

1-A 
b 



e 
d 
c 
b 
a 



At e = this matrix is seven- diagonal. After we set d — it becomes 
pentadiagonal while the next constraint c = makes it tridiagonal. The 
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last, 6 = item in the list of simplifications is the diagonal matrix which is 
compatible with Eq. ffT4|) and non-Dirac, 0^ o )(A) 7^ / at A 7^ 0. 

The most important consequence of our "brute-force" calculations per- 
formed at N = 5 lies in the highly desired clarification of certain higher— iV 
tendencies exhibited by nondiagonal metrics. In the light of this experience 
we were able to formulate and prove the following 

Proposition 4.5. At any dimension N = 3,4, . . . and coupling A G (—1, 1), 
pseudo-Hermiticity relation (TJ7J) is satisfied by Hamiltonian (TJ|) and by the 
bidiagonal pseudometric matrix 



vi N \\) = Q 



(TV) 

(0,(9,0,0,. ..,0) 



(A) 






p 
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p 






(17) 



containing two pairs of X— dependent elements (3 = /3(A) = 1 — A connected 
by two unit diagonals. 



Proof. Equation ffTTj) is compatible with Eq. (J8J 
at N = 4. At any integer N > 4 the insertion of 
into identity. 

Proposition 4.6. At any dimension N = 3, 4 
Hamiltonian (TJJj and pseudometric matrix 



1 at iV = 3 and with Eq. (TT2]) 
ansatz ( fTTj) converts Eq. ( Tl0|) 

□ 

. . . and coupling A G (—1, 1), 
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. 



with two antidiagonals and four elements j3 ■ 
Hermiticity relation (T72P- 

Proof. Equation ( Tl8|) is compatible with Eq. (jSJ 
at iV = 4. At any integer iV > 4 the insertion of 
into identity. 



(1* 



1 — A satisfy the pseudo- 



l at iV = 3 and with Eq. (EE 
ansatz (fl8l) converts Eq. ([10 



□ 



We see that at iV > 3 the same simplification of the interior matrix 
elements to units or zeros occurs in both the bidiagonal and antibidiagonal 
cases. Similar phenomenon will characterize also the structure of all the 



remaining elements of our set of pseudometrics V 



(TV) 



Vf (A). 
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5 The complete set of pseudometrics 



The simplest pseudometric matrices V% (A) entering expansion (Q were 



defined by Eq. (TH|) (where we have to set k — 1), by Eq. (fTTj) (where k = 2), 
by Eq. (JISJ) (where k = N), and by Eq. ([HD (where = N - 1). On 
the basis of inspection of metrics 9^ evaluated at dimensions Af < 5 we 
are now prepared to guess and prove the explicit form of the remaining 
pseudometrics (A) at all the subscripts k such that 3 < k < N — 2. Our 
N hy N candidates c£ for these pseudometrics will be sparse matrices with 
the property 



C 



(N) 



whenever m + n = k (mod 2) . 



(19) 



This means that on a sufficiently large chessboard the nonvanishing matrix 
elements of each of these matrices would only occupy either black or white 
fields. Secondly, all of these "colored" (i.e., black or white) ansatz matrices 
will share the following two-parametric form of their upper segment, 



C 



(AT) 



Cl N \z,v) 






V 


1 






v 


1 


1 



z 


1 



1 







v 


1 


1 






v 



1 





(20) 



copied and shared also by the 7r/2— rotated left segment, lower segment and 
right segment. In compact notation we are led to the left-right as well as 
up-down asymmetric arrays of matrix elements, 



Cl N) (z,v) 





... 

v 

z 

v 

... 







v 



v 




z 

1 



1 


z 





v 



V 










... 

'•• 

v 







z 



. • ' v 

... 
... 



(21) 
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In these two-parametric matrices the uppermost item z sits in the k— th place 
of the first row of our matrix while the bottom line contains z in the (N + 
1 — k)— th place. Together with the other two z's we obtain the quadruplet 
of vertices of connected by the four diagonal rows of elements v forming a 
parallelogram. Keeping still in mind the chessboard visualization of matrices 
Cfc(z,v) we require that this parallelogram separates the outside domain 
(filled just by zeros) from the rhomboidal interior filled by "equal color" 
units and "opposite color" zeros. A few illustrative samples of explicitly 
computed matrices V possessing this structure may be found in Appendix B 
below. We are now prepared to prove our final result. 

Theorem 5.1. Pseudometric matrices V[ N \\) with 3 < k < N — 2 which 
would be compatible with Hamiltonian (TJJ) via condition (TQ) may be identified 
with the X— dependent matrices £j£ (7, 5) where 

7 = 7(A) = ^A * = *(A)~. (22) 

Proof. The inspection of concrete examples using small fixed dimensions 
iV (obtainable by the algorithm outlined in Appendix B) indicates that all 
our elementary pseudometrics Vk{\) may be expected to exhibit the above- 
mentioned fourfold symmetry. This expectation is easily verified by immedi- 
ate insertions confirming that in our algebraic manipulations with Eq. (fit)]) 
it is sufficient to work just with the not too large submatrices ( l20|) glued, if 
necessary, to their rotated neighbors. In this sense, the role of the position 
of the subscript k in interval [3, N — 2] remains inessential. 

In the next preparatory step of our proof let us recall the detailed form of 
our generic ansatz (120]) and visualize it also as written on the black-and-white 
chessboard. Obviously, the color of fields with nonvanishing matrix elements 
will be fixed as "white" or "black" for each subscript k. At the same time, 
our tridiagonal Hamiltonian H will be both "white" (= its main diagonal, 
H (w) ._ 21) an d "black" (= its upper diagonal, H^ b+ \\), as well as its 
lower diagonal, H^ b ~'{X)). This type of coloring simplifies our argumentation 
because in our fundamental Eq. (ITU1) we may use any Pfc(A) = 0^,6) and 
decompose H = H(X) = + H^ b+ \\) + H^(X) and W = H(-X) = 
#H + #(6+)(_A)+#( 6 -)(-A). Obviously, the "color" of any selected C k (7, 5) 
will be shared by its products with and it will differ from its products 
with #( 6± )(±A). Using this idea we may very quickly deduce that starting 
from Eq. (|20|) we shall always obtain product C H in the form characterized 



17 



by its upper segment 
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27 


-7 
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25 -5 
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2 -1 - 



Similarly, the product H * C will be specified by its very similar upper segment 






-5(1 - A) 


27 


-6(1 - A) 





-5 


26 


-1-7(1 + A) 


26 


-5 


25 


-1-5 
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-1-5 


26 


-1-5 
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-2 
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-1 - 
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-2 


2 


-2 


2 



The difference between these two upper-part matrix structures only involves 
the underlined matrix elements. It is easy to check that the coincidence 
of these two remaining triplets of underlined matrix elements is guaranteed 
since it follows from definition fl22|) of quantities 7 = 7(A) and 5 = 5(X). 

The symmetries of our pseudometrics imply that the same coincidence 
of matrix elements will take place for the two lower parts of products C H 
and WC. In contrast, the rotation of the upper-part matrices by the mere 
±7r/2 changes the picture. Different pattern emerges causing, fortunately, 
just an exchange of the underlined matrix elements between the respective 
left or right parts of products C H and C. In this way we arrive at an- 
other, equivalent matrix representation of equation C = C H. Indeed, in 
its new form this representation can be made explicit, say, via its left-part 
embodiment 
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hA) 2 
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We see that the modifications are inessential and that the rotated recipe 
leads to the same identities. We may conclude that the products C and 
C H coincide. This completes the proof of eligibility of our ansatz (j201) + 
for independent non-exceptional sparse-matrix components Vf, (A) of 
the general N— parametric metric at any dimension N. □ 

6 Summary and discussion 

We can summarize that for our square-well model (j5J) the N— term decompo- 
sition (jSj) of all the existing metrics B^(A) may be written in terms of the 
non-exceptional closed-form pseudometrics of Theorem 15.11 complemented 
by the exceptional closed- form metric v[ (A) [defined by Eq. (JHj)] and 
closed-form pseudometrics V% (A) [defined by Eq. (fl7|) ]. V^\{X) [defined 
by Eq. PS])] and Pjf } (A) [defined by Eq. jT5])]. 

In a broader perspective, the existence of remarkable parallels as well as 
differences between the standard and VT — symmetric models of bound states 
has been reconfirmed. It has been emphasized that in the former case one 
works, as a rule, solely with the most elementary Dirac metric Q = I. In the 
majority of VT — symmetric models, on the contrary, the technically most 
difficult task concerns the construction of the appropriate metric or metrics 
8 ^ /. For this reason, the work with the difference (rather than with the 
more current differential) VT— symmetric Schrodinger equations has been 
preferred and may be recommended as easier in technical terms. 

In our present paper we intended to address and formulate the problem of 
solvability in VT— symmetric context. We felt inspired by one of the simplest 
differential square-well Hamiltonians where a "minimal" non-Hermiticity has 
been introduced via boundary conditions and where an exceptionally elemen- 
tary metric has been found in Ref. [TJ. On this background we selected 
and analyzed the difference-equation version (j5J) of this model. 

A supplementary reason for our choice of model (JS]) has been provided by 
Refs. (15] and [2TJ where several physical and dynamical assumptions (say, 
about a large distance between interaction centers, etc) have been made 
in the technically more difficult differential-equation context. In our study 
of the discrete sample of a generic boundary-condition model some of 
the empirical observations made in these references (concerning, e.g., the 
correspondence between the fragile and robust energy levels [15]) reappeared 
and have been illustrated by a few pictures. 

The most interesting conclusions may be extracted from the comparison 
of our present results with their predecessors described in Ref. [T5] . In both 
these cases point-like interactions were used. Still, the decisive advantage of 
their present version has been found in their consequent localization in the 
closest vicinity of the boundaries. Formally, this feature of our H has been 
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reflected by an enormous simplification of the structure of all of the related 
matrices V (= pseudo metrics) and G (= metrics). 

The latter merit of our model was quite unexpected and its explanation 
also forms a mathematical core of our present message. A posteriori we 
may conclude that the simplification of our matrix equations for V and 
resulted from a "hidden" possibility of their split in the "white" and "black" 
components. This is also a deeper reason why the simplicity of the metrics 
attached to our present model is in a sharp contrast with the complicated 
recurrent nature of the analogous matrices in the older model of Ref. |13j . 

In connection with this feature of our model the core of feasibility as well 
as of the rigorous form of our constructions may be seen in the availability 
of appropriate ansatzs. They were found by extrapolation from investiga- 
tive constructions performed, at the smallest dimensions, by the brute-force 
linear-algebraic techniques. A posteriori we must appreciate, therefore, the 
drastic reduction of the large set of algebraic quasihermiticity conditions (J7J) 
or ([TO]) to the mere double definition (122]) of functions 7(A) and 5(A). This 
reduction seems to have been caused by a certain purely formal interplay 
between the tridiagonality of H and its free-motion character preserved near 
the origin. It was precisely this fine-tuned dynamical input which suppressed 
the computational difficulties and which facilitated, decisively, the explicit 
interactive and extrapolative analysis of Q(H). 

Our closed formulae appear transparent, especially if one decides to work, 
say, with just a few terms in the general series ([6]). In this way even the very 
pragmatic users of the discretized norm or inner product 

«V#> := (m (S) = EE r(xk)®k,n<P(x n ) in 

k n 

might employ not only the diagonal matrix metric of Eq. (TH|) but also, 
say, a Sobolev-space resembling discretized inner product with tridiagonal 
q(n) _ 2v[ N ^ — ^V^" 1 where, even at large N, the obligatory guarantee of 
positive definiteness would just require that I7I < 1. 

In a more general framework of study of mutually non-equivalent possible 
physical Hilbert spaces assigned to a given Hamiltonian H we believe 
that a deeper insight if not classification could be obtained in the nearest 
future, especially on the ambitious level aiming at the less elementary inter- 
actions and/or more complicated combinations of independent components 
in the positive definite matrices of metrics. 
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Appendix A: VT— symmetric models 



One of the most common tacit assumptions which lies at the very heart of the 
correct physical interpretation of Eqs. and/or (J2]) is that the real variables 
x and/or Xk represent an experimentally measurable quantity (most often, 
a coordinate or momentum of a point particle). This type of postulate has 
been declared redundant in the so called VT — symmetric models where all 
x or Xk are allowed complex (cf. Refs. [3] for a full account of the theory). 
The endpoints ±L may be then chosen as any left-right symmetric pair of 
points, finite or infinite, in complex plane of x. 

This is a new freedom. It implies, e.g., that the "wrong-sign" potential 
V(x) = —x 4 becomes tractable as a fully legal source of a discrete and real 
spectrum of bound-state energies which is bounded from below. One only 
has to complexify the points of boundary in an appropriate left-right sym- 
metric (called, for historical reasons, VT — symmetric) manner. In a typical 
replacement +L — > gexp(-iip) and — L — > — gexp(+i(p) one uses a very large 
(or infinite) real g ^> 1 and some safely nonvanishing real angle <p (cf., e.g., 
Refs. [22] for more details). 

In the VT — symmetric scenario the role of the variable x (and, mutatis 
mutandis, of Xk) is purely auxiliary. In principle, this variable does not 
represent an eigenvalue of any operator of observable even when it remains 
real. Mathematically one speaks about a "false" Hilbert space equipped 
with the most common definition of the inner product of wave functions, 

(Val^) HV>«I^) (F) ) = / ttWMtfdx in H {F) . (23) 

Hamiltonians are represented there, typically, by non-self-adjoint differential 
or difference operators, 

H = — f- + V(x) + # f in HS F) . (24) 
oar 

As long as they have to generate a unitary time evolution, we must change 
the definition of Hermitian conjugation in order to make them properly Her- 
mitian in the resulting "standard", i.e., physical Hilbert space 

H = —^ + V(x) = H t (= e-^e) in H {s) . (25) 

The auxiliary operator O = O^ > responsible for such a change of the 
definition of Hermitian conjugation in Ti.^ is called "metric" |1Q| . All the 
other operators O of observables must be "cryptohermitian" , i.e., Hermitian 
in the same space, 

o = o x (= e^o+e) in n {s) . (26) 
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The latter requirement (or a set of requirements if necessary) specifies the 
physics of the system in question. In opposite direction, the set of all the 
"hidden hermiticity" requirements (|25p + (T2"oT) may be understood as a prac- 
tical recipe for the explicit determination of the correct and unique metric 
operator 6 assigned to given H and O [TO], [17] . 

An immediate consequence of the latter formulation of the theory is that 
the two Hilbert spaces H {F) and H {s) coincide as vector spaces composed 
of the same wavefunction elements. The main and only difference between 
them lies in the update of the inner product in the latter space. In Ref. p2] 
a compactified notation using a double-bra symbol has been recommended, 

(Va|V>&> =^ (V>a|e|V&> (= m^) iS) = ((V'ai^)) • (27) 

Our Hamiltonian H and Schrodinger Eq. (or, mutatis mutandis, Eq. (T5]), 
plus all the other operators O of observables) find their physical probabilistic 
interpretation in the Hilbert space Ti^ with a nontrivial metric B. Thus, 
the internally consistent concept of the solvability of the model should, in 
principle, involve not only the feasible construction of all the wave functions 
ip(x) and of all the related bound-state spectrum of energies E but also the 
practical feasibility of the assignment of the metric G to our system. This is 
a challenging problem, addressed in our present paper. 



Appendix B: The sample of computation of 
all the pseudometrics in the six-dimensional 
square-well model (0) using the linear set of 
Eqs. (HUD 

The twelve unknown parameters entering the symmetry-reduced N = 6 
ansatz for the pseudometric 

a b c d e j 
b f g h k e 
c g m n h d 
d h n m g c 
e k h g f b 
j e d c b a 

must be shown compatible with the thirty six quasi- Hermiticity requirements 
represented by Eq. (TlOl) . The reduced set of six items of these equations (say, 
number 2, 9, 3, 10 and 4) offers the affirmative answer. Indeed, this sextuplet 
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of equations 

-f + f X + a + aX + c = 
— c — c A — m + f + h = 

-g + gX + b + d = 
—n — d — dX+g+k=0 

-h + hX + c + e = 
-k+kX+d+j+jX=0 

leaves six parameters unconstrained and, with abbreviations 
u = c + a(l + X) + cj^ _ c + = a (1 + A) + c (1 + A 2 ) 



A 



1-A 



1-A 



1-A 



1-A 



and 

^ b + d d + j(l + X) J/1 x . j(l + A) d(l + A 2 ) 
it specifies the following complete solution of Eq. (flOl) at AT = 6, 



1-A 



(6) 

(a,6,c,d,e,i) 



a 
b 

c 
d 
e 
J 



b 

c+a(l+X) 



1-A 

b+d 

1-A 
c+e 

1-A 
d+j (1+A) 



C 
b+d 
1-A 

17 



d 

c+e 
1-A 

7 



e+e y b+d 



1-A 

e 



c+e 
1-A 

d 



b+d 
1-A 

C 



1-A 

c+e 

1-A 
b+d 
1-A 
c+a(l+A) 

1-A 

b 



J 

e 
d 
c 
b 



At j = this matrix is nine- diagonal. It further becomes seven-diagonal 
with e = 0, pentadiagonal after fixing d = while the next constraint c = 
makes it tridiagonal. Finally, at b = we arrive at the simplest, diagonal 
metric predicted by Eq. (I14p . Due to the symmetries of our problem we may 
infer that, in parallel, there also exists a very similar antidiagonal metric (II 5p 
containing, again, all units up to the endpoint exceptions. 

The similar simplifications of the "interior" matrix elements in V^(X) 
may be noticed to appear in the tridiagonal and antitridiagonal cases where 
the two particularly elementary special pseudometrics given by the respective 
Eqs. f|T71) and ffTSl) with (3 — 1 — A are reproduced, 
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The remaining two independent components of the set of pseudometrics come 
out as parametrized via 7 = (1 — A)/(l + A 2 ) and 5 = 1/(1 + A 2 ) yielding 



n b> w = 



7 

s s 

7 1 5 

5 1 7 

5 5 

7 



7 

5 5 

5 1 7 

7 1 5 

5 5 

7 



in full compatibility with Theorem 15.11 
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